Multiple Quantum Phase Transitions of Plutonium compounds 
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We show by quantum Monte Carlo simulations of realistic Kondo lattice models derived from 
electronic-structure calculations that multiple quantum critical points can be realized in Plutonium- 
based materials. We place representative systems including PuCoGas on a realistic Doniach phase 
diagram and identify the regions where the magnetically mediated superconductivity could occur. 
Solution of an inverse problem to restore the quasiparticle renormalization factor for /-electrons is 
shown to be sufficiently good to predict the trends among Sommerfeld coefficients and magnetism. 
Suggestion on the possible experimental verification for this scenario is given for PuAs. 

PACS numbers: 71.27.+a, 75.30.Kz, 75.40.Mg 



Motivation Discovery of unconventional supercon- 
ductivity in PuCoGas [l|, Q opened a new arena for the 
studies of strongly-correlated materials. It has the high- 
est superconducting transition temperature T c = 18.5 [K] 
among /-electron-based materials and it has been dis- 
cussed to reside somewhere in between the Cerium- 
based heavy fermion (HF) superconductors and high-T c 
cuprates, where the latter still challenges theoretical con- 
trol from first-principles. 
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In the present work we make predictions on the 
magnetism and HF behavior of several Pu compounds 
including Pu-115's where the mechanism of possible 
magnetically-mediated superconductivity is discussed to 
be more complicated than their Cerium counterparts Q ■ 
Furthermore, with experimental challenges such as the 
self-heating of samples due to the radiative nature of Pu 
nuclei, a computational guide should be of help regarding 
the determination of the linear coefficient of electronic 
heat capacity, so-called Sommerfeld coefficient 7. Our 
computational method is based on a recently-developed 
scheme for realistic Kondo lattice simulations 

as 

which 

enabled us to predict the location of magnetic quantum 
critical point (QCP) 6] from electronic structure calcu- 
lations 7} for HF materials. 

Our main results arc shown in Fig. [T]in the format of 
Doniach phase diagram Q plotted with realistic settings 
for the target materials. Striking double-dome structure 
is seen both in the magnetic Doniach phase diagram of 
PuCoGa 5 plotted on ( Jk, T^)-p\a,ne and an analogous 
plot on ( Jk, y)-plane, where Jk is the Kondo coupling 
and Tn is the Neel temperature. For Pu-115's there are 
at least two antiferromagnetic phases with at least three 
QCP's. We see that Pu-115's are indeed separated from 
the first i.e. lowest-energy magnetic QCP, being consis- 
tent with the situation discussed in Ref. |3| . However we 
find that the second and third QCP's are encountered on 
JK-axis and the realistic point for Pu-115 is actually in 
the proximity to the latter. 
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FIG. 1: (Color online) Summary of our results for Pu-115's 
on realistic Doniach phase diagram. The realistic data point 
is indicated by the symbol on the line for each target ma- 
terial. The main panel show the trend among Sommerfeld 
coefficients and the inset show the Neel temperatures. 



Methods The electronic structure calculation based 
on local density approximation (LDA) combined with dy- 
namical mean-field theory (DMFT) (LDA+DMFT) has 
been successful in addressing interesting properties of 
strongly-correlated materials [7|. What motivates us for 
the Kondo-lattice model (KLM) description of HF ma- 
terials is that efficient and exact quantum Monte Carlo 
(QMC) simulations in low-temperature region are possi- 
ble Q, typically around O(10) [K] and down to 0(1) [K]. 
This advantage is due to having only /-spins and elim- 
inating the /-electron charge deg rees of freedom via 
SchriefFcr- Wolff transformation [10( implemented in a re- 
alistic way This is in contrast to the fact that the 
standard LDA+DMFT based on solving the Anderson 

i3 for <5-Pu, 



model, which was used e.g. in Ref. [11|J for d-Pu, typi- 
cally can reach the temperatures down to O(100) [K] if 
the core impurity problem is to be exactly solved by QMC 
method. Here, some basis-cutoff schemes have been im- 
plemented to reduce the computational cost. 

One of the reasons Plutonium compounds have been 
interesting and difficult to address is that they reside on 
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material iV tot (0) jV-^O) -TrSA(0)/7r [eV] 



PuCoGa 5 


34.33. . 


1.203. . . 


0.705 


PuRhGa 5 33.03. . 


1.494. . . 


0.912 


5-Pu 


20.40. . 


2.781... 


1.13 


PuSe 


17.30. . 


4.342. . . 


0.619 


PuTc 


13.45. . 


0.7729. . . 


0.360 


PuAs 


8.351. . 


0.6068. . . 


0.255 


PuSb 


10.01. . 


0.2699. . . 


0.193 


PuBi 


7.723. . 


0.2534. . . 


0.133 



TABLE I: Summary of LDA+Hubbard-I results for target 
materials. The unit of density of states, N tot (0) for all 
electrons and N^(0) for /-electrons on the Fermi level, is 
[states/Ry/cell]. 



the border of itinerancy and localization of 5/-electrons 
among actinides III, 14 1. At least for Pu-chalcogenides 
and pnictides, experimental evidence for localized 5/- 
electrons was revealed [H] and a recent theoretical 
work [lf| agrees with that so these can be benchmark 
cases for the realistic KLM simulations. For Pu-115's 
it has been known that the Curie law persists down to 
the superconducting temperature for PuCoGas [l[ which 
supports the presence of localized 5/-electrons, but some 
attention must be reserved for a possible sample depen- 
dence: radiative Pu decays into U, which can introduce 
magnetic impurities. It is thus controversial whether the 
Curie- Weiss law is intrinsic or not Analyses of ex- 
periments point to n/ = 5.03 for <5-Pu [13] which we 
believe is sufficient for the KLM to work. The valence de- 
duced from the calculations shows a much larger spread 
(between 4 (3 and 6 0) with the values of 5.2 for 6- 
Pu [2(| and 5.26 for PuCoGa5 [2l[ from most accurate 
CT-QMC calculations. 

Our realistic KLM framework for the above-mentioned 
Pu compounds goes as follows. At the first stage LDA for 
s, p, and d-conduction electrons and Hubbard-I approx- 
imation 22| for the self-energy of localized /-electrons 
gives us the partial densities of states and hybridiza- 
tion functions as prescribed by LDA+DMFT frame- 
work [7]. The data are summarized in Table HI It 
is clear that Pu-115's have much higher energy scales 
than the Cerium ones At the second stage we 

solve the low-energy effective KLM Hamiltonian with 
dynamical-mean field theory 



23l . |24| |. utilizing state- 



of-the-art continuous-time quantum Monte Carlo (CT- 
QMC) impurity solver 0, - For the 5/-orbitals of 
Pu, it is known that there is a big spin-orbit splitting 
of 1 [eV] and the five possibly localized electrons fill in 
the lower j = 5/2 multiplet up to leaving one localized 
hole [2(| . Photoemission experiments show that the level 
of the localized hole is separated from the Fermi level by 
1 [eV] which is verified by our theoretical estimates. We 
neglect the crystal-field splittings which are known to be 
small in Pu compounds in the local 5/ level. 



Solving an inverse problem to restore /-electrons 
Even if we eliminated /-electrons and kept only /-spins 
in our KLM, part of the information for the localized 
/-electrons can be restored from the relation E c (iu; n ) = 
V 2 /[iu! n — e/ — £/(mj„)], where E c is our conduction- 
electron self energy, ioj n — (2n + l)irT is the Matsubara 
frequency, e/ = — 1 [eV] is the position of the impurity 
level, and E/ is the /-electron self-energy which we do not 
have explicitly in our KLM calculations. Provided that 
we reach the temperature for a given target material to be 
in a Fermi-liquid region concerning its /-electrons, which 
is mostly the case for Pu compounds, the quasiparticle 
renormalization factors are well defined and written as 
z x = (1 — d^sY, x (ioj n ) j d(iuj n )) 1 , with x — c and / for 
conduction electrons and /-electrons, respectively. Here 
the derivative is taken at iui n = 0. We get from the 
above definition of S c the following inversion relation : 
Zf = [|S c (0)| 2 /y 2 ]z c /(l — z c ). Because z's are written in 
terms of the derivative of the corresponding self-energy 
at the lowest frequency, our effective low-energy descrip- 
tion based on KLM enables a good solution of this in- 
verse problem as far as Zf is concerned. The Sommerfeld 
coefficient 7 = (l/3)7r 2 A' r ff (0), where N c g(0) is the effec- 
tive total density of states (DOS) on the Fermi level, can 
be estimated by N cS (0) = N c (Q)/z c + N f (0)/z f , where 
N c (0) is DOS of s, p, d-conduction electrons and Nf(0) 
is that of localized /-electrons in our LDA+Hubbard-I 
calculations. With a given KLM, we extract z c from S c 
obtained after DMFT, invert it to Zf, and get the Som- 
merfeld coefficient 7 with the above formula. In this way 
we can restore an analogue of Doniach phase diagram for 
7 as was shown in Fig. Q] for Pu-115's. The results on 
the realistic data point for each target material are sum- 
marized in Table [Til together with the experimental data 
taken from the literature. Our prediction follows the ex- 
perimental trend among 7 semi-quantitatively. We note 
that 7 is sensitive to the estimate of the realistic point of 
Jk especially around QCP's, considering the sharp peak 
structure as seen in Fig. [TJ for the plot of 7 vs Jk- So 
the overall trend among materials is the most important 
result. 

Magnetism and quasiparticle renormalizations The 
results for magnetism are schematically summarized in 
Fig. [2] for all target materials in the format of a rescaled 
Doniach phase diagram. It illustrates how we understand 
the results in the inset of Fig. [TJ for Pu-115's. Striking 
multi-dome structure shows up together with multiple 
QCP's for materials with strong Kondo coupling. We find 
that Pu-115's are located in a region where antiferromag- 
netic long-range order is suppressed, possibly near a hid- 
den or pseudo-QCP, within some numerical noise at the 
lowest reachable temperatures at present. Inspecting the 
distribution of materials around the QCP's in Fig. [2j we 
have pnictides on the left-hand side and chalcogenides on 
the right-hand side of the antiferromagnetic QCP. This 
is consistent with what has been known experimentally, 
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TABLE II: Summary of our data obtained with realistic 
Kondo lattice simulations and our prediction for 7 based on 
them. The unit of 7 is [mJ/mol/K 2 ]. Experimentally known 
results are taken from the literature. 
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FIG. 2: (Color online) Schematic summary of our magnetic 
phase diagrams for Pu compounds plotted on the (t, TN)-plane 
where t = (Jk — Jk,qcp)/Jk,qcp is the rescaled Kondo cou- 
pling with Jk.qcp being the first QCP in (a) and the third 
QCP in (b). 



that is, pnictides such as PuAs, PuSb 0, and PuBi [HI 
arc magnets and chalcogenides such as PuSe and PuTe 
are paramagnets [2^]. The actual magnetism is strongly 
spatially anisotropic 3 IH whose treatment is unfortu- 
nately beyond the level of single-site DMFT description. 
For now we will leave the issue of ordering wavevectors 
for future projects and focus on the trends across target 
materials spanning between magnetism and HF behav- 
ior. The characteristic energy scales of Kondo-screening 
and magnetic ordering have been captured by fully incor- 
porating the frequency-dependence of the hybridization. 

Multi-dome structure together with multiple QCP's 
shown in Fig. [T]for Pu-115's can be understood in terms 
of strong-coupling nature of the Kondo lattice 36j based 
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3: (Color online) Analogue of Doniach phase diagram 
a) Kondo temperature and (b) quasiparticle renormal- 



FIG 

for ( , 

ization factors for PuCoGas. The inset in (b) is a zoom-up 
picture around the first QCP with the vertical axis plotted in 
logarithm. 



on the growth of characteristic Kondo energy scale Tk 
with respect to Jk as shown in Fig. [3] (a) which is ob- 
tained from our local susceptibility data. Simple pertur- 
bative arguments also prompt for two crossover points 
between Trkky ~ JkP ano - -^K ~ ex P[ — 1/(^kp)]/p tem- 
perature scales, where the latter can saturate at some 
point with respect to large Jk while the former keeps on 
growing. Here p is the characteristic DOS. 

We demonstrate our predictive power regarding Xk 
also for the case of 6-Pu where we get 7k ~ 10 3 [K] from 
our local susceptibility data which is seen to be close 
to the previous results, T K - 700 [K] in Ref. 0. We 
note that the Kondo screening energy scale is approach- 
ing a comparative scale to the characteristic bandwidth, 
or the kinetic energy of the conduction electrons which 
is O(l) [eV 1 . Such situation had been discussed in the 
literature [3Q . in the context of models, which is now 
found to be realized in Plutonium heavy-fermion materi- 
als. 

The behavior of quasiparticle renormalization factors 
z x (x = c or /) as shown in Fig. [3] (b) further gives 
the physical picture: starting from Jk — where there 
are free conduction electrons and completely localized f- 
electrons with (z c ,Zf) = (1, 0), the former gradually gets 
renormalized toward the first QCP and /-electrons gets 
"delocalized" in the sense that they start to take part 
in Fermi surface (FS) [H, Passing the first QCP, 
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heavy quasiparticles composed both of conduction elec- 
trons and /-electrons evolve together after z c has shown a 
dip around the first QCP, letting /-spins show up again 
with the underscreening effects that correspond to the 
slightly elevated z c . Thus after the revival of magnetism 
the same thing can happen again and could repeat it- 
self, with the re-defined much smaller energy scales every 
time, all the way to the Jk ~^ oo limit, being consistent 
with the phase diagram obtained in Ref. |36j |. 

This strong-coupling KLM scenario for Pu compounds 
can in principle be checked by de Haas-van Alphen ex- 
periments which would measure the size of the FS to 
see if it counts the number of "localized" /-electrons. 
The ferromagnetic phase of PuAs should be the one 
with the "large" FS including the spins of localized 5/- 
electrons. This phase would be in contrast to the typi- 
cal magnetic phases in HF compounds with the "small" 
FS, being located in the weak-coupling region. Following 
the method of Ref. [3^, we can track the evolution of 
large FS for representative Pu compounds obtained from 



-3?£ c (jw„) 



iui n =0 



and we find that PuAs shows a remark- 



able evolution of "large" FS, which is to be compared 
with experiments to see if the strong-coupling KLM pic- 
ture can hold. 

Conclusions We have found that multiple QCP's 
show up for Plutonium-based materials that have 
stronger Kondo couplings than their Cerium counter- 
parts. Our methodology captures the quasiparticle renor- 
malization factor and the characteristic energy scale cor- 
rectly. The striking multi-dome feature of the Doniach 
phase diagrams for Pu-115's and 5-Pu as well as the mag- 
netic and HF behavior among Pu pnictides and chalco- 
genides is understood on the basis of strong-coupling 
limit of KLM. This picture can be verified by looking 
at the size of the Fermi surface for PuAs. 
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